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Abstract 

Spectrum of gluons in the adjoint source field is computed ana- 
lytically using the QCD string Hamiltonian, containing only one pa- 
rameter - string tension, fixed by meson and glueball spectrum. Spin 
splitting is shown to be small. A good agreement is observed with spa- 
cially generated gluelump states measured on the lattice. Important 
role of gluelumps defining the behaviour of field-strength correlators 
is stressed and correspondence with earlier computations of the latter 
is established. 



1 Introduction 

Gluelumps are not physical objects and their spectrum cannot be mea- 
sured on experiment. However, as will be discussed below, they play more 
fundamental role than any other hadrons since gluelump masses define field 
correlators in the QCD vacuum and in particular string tension. The purpose 
of the present study is to calculate analytically the gluelamp spectrum in the 
framework of the so-called QCD string model (QCDSM). This model was 
developed in the earlier papers [|, f| for spinless quarks, and was augmented 
by the spin-dependent terms treated as a perturbation 0, [J . The resulting 
Hamiltonian of the rotating string with quarks or gluons at the ends was 
recently systematically applied to mesons, hybrids and glueballs (see review 
in H and recent developments in |7|, §). Three features are characteristic 
for the model. First, QCDSM is directly derived from QCD with few as- 
sumptions supported by the lattice data, i.e. the minimal area law for the 
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Wilson loop, and the dominant role of valence quarks and gluons. Second, 
it is fully relativistic and the Hamiltonian can be obtained in the cm. sys- 
tem m, or in the light-cone system || or else on any other hypersurface. 
Third, the model contains the minimal number of parameters: fundamental 
string tension 0/ = 0.18 GeV 2 defined by the meson Regge slope, the strong 
coupling constant a s (taken as constant in first approximation, a s (/i), where 



\x is an inverse size of the system, or with the freezing behaviour a s (r) [10 
in more accurate calculations), and finally the quark (antiquark) selfenergy 
Co = —0.25 GeV, which should be subtracted for each quark from the mass 
of a meson, hybrid, or baryon. Note, that for gluons this term is absent be- 
cause of gauge invariance and therefore spectrum of glueballs and gluelumps 
is defined only by Of and a s . 

It is remarkable, that with these fixed universal parameters one obtains 
spectrum of all hadronic systems in good agreement with lattice data and 
experiment 0. In particular, the glueball spectrum depending only on Of 



(with spin splittings depending also on a s ) comes out [ I I ] in remarkable 



agreement with recent lattice data [12 



At the same time the heavy-light mesons have been calculated in this 
method in [[L3|] including quark decay constants /m, M = B, B s , D, D s again 
in good agreement with experiment and other approaches. 

This gives a hope that the QCDSM can be used also for the system 
containing a valence gluon connected by an adjoint string to the adjoint 
source ("infinitely heavy gluon") - the so-called gluelump [T|, EJ. 

Gluelump has a resemblance both to glueballs and to heavy-light mesons. 
In the framework of the string model, gluelump is a glueball with the adjoint 
string between gluons when one gluon is made static, and one expects the 
same structure of the spectrum of the LS-coupling type, which is obtained 
IT I and similar to the lattice spectrum calculated in fH2f . 



111 



On the other hand, gluelump differs (within the model) from heavy- 
light mesons only in that the fundamental string is replaced by the adjoint 
one (plus minor differences in the spin splitting terms and the role of color 
Coulomb forces). 

Therefore the calculation of the gluelump spectrum is expected to be as 
successful as in the case of glueballs and heavy-light mesons. From physical 
point of view the gluelumps are important at least in two respects. First of 
all gluelumps play a fundamental role in the nonperturbative structure of the 
QCD vacuum, since gluelump Green's functions are connected to the field 
correlators (F flu (x)Fx a {0)) or its invariant parts D(x), Di(x), introduced in 
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n|, which are basic elements in the Method of Field Correlators (MFC). 
In particular, the gluon correlation length T g , defining the nonperturbative 
dynamics of confinement (I5| , is simply the inverse mass of the lowest 
gluelump, and was computed in this way before in several systematic studies 

Secondly, gluelump masses define the screening length for the static po- 
tential in given representation D of SU(3), namely the minimal length of the 
string, which may decay into two gluelumps. 



Recent precise measurements of static potentials Voir) made in (18 



vealed an accurate (within 1%) Casimir scaling without any signs of screening 
(string breaking) at all distances 0.05 fm < r < 1.2 fm. To understand these 
results one should know screening length for all D. 

The structure of the paper is as follows. In the first part of the paper 
the gluelump spectrum is computed, with the corresponding Hamiltonian 
constructed in section 2, spin splittings considered in section 3, and spec- 
trum calculation in section 4. Comparison to the lattice data and the bag 
model is given in section 5. Section 6 contains conclusions and outlook, 
while Appendix is devoted to the details of hyperspherical formalism used 
for calculation of two-gluon gluelumps. 



2 The gluelump Hamiltonian 

The starting point is the gluelump Green's function G glump (x, y) which is 
obtained from the initial and final gluelump operators \]>( m > ou *) ) expressed in 
terms of the valence gluon field and background gluon field so that the 
total field is = B^ + a^, with the gauge transformation properties 

B^U+(B^ + -d^)U, a,^U + a,U. (1) 

Similarly to the glueball case one can write 

G 9ium P ^ y) = {t r adj (T(° ut \x)G^(x,y)T m (y)$(y,x))) B , (2) 

where Q(x,y) is the parallel transporter represent at ing the adjoint source 
Green's function: 

$(x,y) = Pexj>(ig [* B^dzJ, B^ = B^ (3) 
Jy 
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and G^ u (x,y) is the valence gluon Green's function || 

G„ u (x,y) = (D\ u + 2igF ia ,(B))-* = 

= j o ds J Dze- K Pexp[tg J B fl dz fl + 2g ^ F^(^(r))rfr] (4) 

In (Q), r( ott *' m )(x) are operators defining quantum numbers of final and initial 
gluelump states, below in Table 1 listed are examples for the lowest mass 
states. 

The last term in the exponent in (ffl) describes the interaction of the 
gluon spin with the background (and with the adjoint string, which will 
be created out of it). Here as well as in the glueball case |TTJ we shall treat 
this interaction as perturbation to be discussed in section 3, and disregard 
it in the first approximation. Then Eq. (g) reduces to the average of the 
adjoint Wilson loop with the contour C consisting of the straight line of the 
source and the path z{r) of the valence gluon, integrated upon in (f|). 

Assuming the area law for the adjoint Wilson loop in agreement with 
lattice measurements [Q, one can go over in a standard way P, |4| from 
the Green's function G^ 9lump ' to the corresponding Hamiltonian H <y9lump \ 
connected by 

G 9lump (x,y) = (out\exp(-H {9lump) \x -y\)\in). (5) 

The resulting Hamiltonian can be obtained from that of the qq system 0, || 
when mass of one quark is zero and mass of another going to infinity, while 
the quark string tension aj is replaced by a ad j = |cr/. One obtains 

H (glump) = ^ + rf + L ( L + l )l r2 

"2 2// 2( Ji + fidP0*v(J3)) 

Here \i and v{0) are respectively the constituent gluon mass and energy 
density along the string coordinate /3, to be found from minimization of the 
Hamiltonian. 

For the angular momentum L = the stationary point in v(/3) yields a 
simple answer 

H<*™p\L = 0) = £ + £ + * a4i r. (7) 
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For L > instead of solving for the complicated operator (Q) we shall treat 
effect of L perturbatively as in ||], yielding accuracy around 5% for L < 3, 
namely 

2 

H ( 9 ium P )( L >0) = f + | + a a4?r + AM L (8) 
where we have defined as in [|] 

and /i in @ is to be found from the minimum of the mass eigenvalues of 
fj(giump)^ ( accuraC y f this procedure was found in || 19 1 to be better than 
5%). 



3 Spin-dependent and color Coulomb inter- 
action 

We now take into account the term F^ u in which yields the gluon spin S 
dependence, since 

- 2iF ik = 2SB, ik = 1, 2, 3. (10) 

The analysis of the spin-dependent terms can be done as in the case of 
glueballs ]TTJ] and heavy quarkonia || (note however, that nowhere we use 



the nonrelativistic inverse mass expansion, instead the lowest (Gaussian) field 
correlator (FF) is retained in the average of (f|), this procedure was tested 
20| to yield accuracy around 1%). 



in 



As a result one obtains 

AH LS = AH LS {Thamas) + AH LS (pert) (11) 

where 

AH LS {Thomas) = -^—(^ = ) LS (12) 

and f(r = oo) = 1 takes into account that confining interaction starts at 
small r < T g , as r 2 , which significantly decreases the matrix element {^-), 



see for details, and |Tl[ for numerical results for the glueball case. 

With AH(pert) the situation is more subtle. Indeed the one-gluon ex- 
change does not contribute to the potential V\(r) (we stick to the standard 
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notations of Eichten-Feinberg-Gromes approach), and the latter is the only 
LS piece surviving for the gluelump case. Hence one should consider the 
one-loop contribution, similar to the one considered in [21] for heavy quarko- 
nia, in the limit of one mass infinite and another equal to fi. In absence 
of actual calculations for the gluelump case, and to have an orientation, we 
shall translate the result of [[H] replacing Cp — > Ca, and mi — > fi, m 2 — > oo. 
The result is 

T S ry 2 C 2 

AH LS (pert) = 2 * A (2 -hi fir- J E ) = 6 (per<) LS. (13) 

As we shall see the overall magnitude of AHls and of its parts is small, 
and finally can be neglected. For an estimate we shall take the brackets in 
(P~5| ) to be equal to one, and the matrix element (^) can be found from the 
equality [|TT1 

L(L+1)(I) =fia adj (14) 

Consider now the contribution of perturbative gluon exchanges to the inter- 
action between the valence gluon and static source. 

Here we adopt the same attitude that was suggested and discussed in 
detail in |l]J . Namely, the perturbative ladder for the gluelump is equivalent 



to the BFKL ladder of the pomeron (more exactly to one-half of that ladder 



cut along the valence gluon line). For the BFKL ladder it was shown [22 
that the higher-order correction practically cancels the lowest order result. 
This fact was used in the glueball case in [T| to neglect the Coulomb ladder 



completely, and the resulting glueball masses agree surprisingly well with 
the lattice data (whereas retaining color Coulomb interaction would reduce 
the lowest masses by 0.5 GeV, and increase spin splittings 3 times in strong 
disagreement with lattice data). 

Accordingly for the gluelump case we at first also disregard color Coulomb 
exchanges both in mass eigenvalues and in wave functions, so that relation 
([14]) holds true and the matrix element (^p-) is to be calculated with eigen- 
functions of (|7|) and (^). One finds that the spacial size of gluelump (in 
absence of color Coulomb force) is y/2 times less than that of two-gluon glue- 
balls, and is around 0.5 fm for lowest states. Hence the influence of (however 
reduced) Coulomb force in this system can be important. 

Therefore also calculations have been done for two other values of effective 



a s , which when taking into account higher loop effects [22] reduce to a smaller 



value, denoted by a s . Two values of a s , a s = 0.15 and a s = 0.195 have 
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been used and the results for the spin-overaged mass are shown in Table 2, 
while the values of spin-orbit matrix elements h^ Thomas ^ and &( peri ) change for 
a s = 0.15 at most 1.5 times and not listed in the Table. 

The resulting values of spin-splitted masses are given in Table 4 for a s = 0. 



4 The gluelump spectrum 

We are now in position to calculate the spectrum of (H) and the spin splittings 



The total mass eigenvalue is written as 

M(J,L,n r ) =M (n,L)+AM L + AM LS (15) 

where M (n, L) is the eigenvalue of the operator h(fjL) = fJ,/2 + p 2 /2/j l + cr ad jr, 
minimized over the values of \x. 

h(^(x)=e(fi)^(x). (16) 

Using the standard technic of the QCDSM || ||, one finds 

=2 + MVS 

where a(n, L) are the eigenvalues of the reduced equation found numerically 
and given in Table 1 (lower entries). Minimization of ( |T7D over fi yields for 
/i = (jLq the values given in Table 1 (upper entries). One has from ([FT 



M^) = (%^) 3/4 (2^) 1/2 . (18) 

From e(n) (|17]) one immediately finds that 

M (n, L) = 2^ (n, L) (19) 

The values of Mo(n, L), AMi and the resulting spin-averaged masses M(n, L) i 
Mo(n, L) + AMl are given in Table 2. 

Now we turn to the spin splittings of the obtained levels. The values of 
(■\) are taken from (14) and those of (^-) from the corresponding glueball 
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matrix elements in Table 7 of ||llj| , with the proper rescaling from the two- 
gluon Hamiltonian to h(n) The resulting splittings AM^s can be written 

as 

AM LS = J( J + 1 )- L ( L + 1 )~ 2 ^(Thomas) ^ ^ + ^ert) ^ L)) (2Q) 

The values of U Thomas ^ and U pert ^ are given in Table 3. Finally summing up 
all corrections one obtains the resulting gluelump masses listed in Table 4 
for a s = 0. 

Now we turn to the two-gluon gluelumps, which correspond to C = +1. 
This is a special case of 3g glueballs when the mass of one gluon tends to 
infinity. Therefore one can use the same Hamiltonian technic as was exploited 
for 3g glueballs in |IT|, i.e. the Hamiltonian h(fi) in ( fIB| ) is replaced by 
h(fii, fiz) (with /ij to be found again from minimization of eigenvalues). 

hfato) = y + y + ^ + ^ + Of + |r 2 | + |r x - r 2 |} . (21) 

Note the appearence of <jf in (0), since the string is now fundamental 
and forms a triangle, passing through two gluons and the adjoin source. 
The Hamiltonian ( pl|) can be treated using hyperspherical formalism, for 
details see Appendix. Here we only quote the result for spin-averaged masses 
(obtained with accuracy better than 5%) 

1 K(K + 4) 

M(n, K) = 4/i + u (n + -) + 1 ^ M /i (22) 
where n — 0, 1, 2, ...;K = 0, 1, 2, ... and 

/i = (J^ ' (V2 + 2) 1 / 2 , Wo = -^/i . (23) 

For the lowest state, K = 0, n = one obtains 

M (29) (0, 0) = 2MGeV(a f = 0.18GeV 2 ). (24) 

The hyperfine splitting of this level into J PC = ++ and 2 ++ can be calculated 
in the same way as in [|ll| with the result that ++ lies below 2 ++ with 
the interval proportional to a s (see Appendix for explicit expression). We 
place the spin-averaged value (|24]) in Table 4 and hyperfine-splitted level 
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++ in Table 5. It is interesting to note that in two-gluon gluelumps both 
hyperfine and tensor forces are operating in contrast to one-gluon gluelumps, 
and therefore the mixing of the state B^B*- 2 - 1 (which is orbitally excited 
with respect to a'paf'^) with the latter is significant, which may explain the 
appearence of this state in lattice calculations. 

5 Discussion of the spectrum. Comparison 
with lattice data 

The spectrum given in Tables 2 and 4 demonstrates the following features. 
First of all one can see that the spin splittings of levels are small and can 
be neglected in first approximation. The ordering of levels has the same 
character as for two-gluon glueball masses and can be roughly described by 
the equation 

M(n, L) = (2n + L)u + M (25) 

where M (a s ) (1.5 - 2.5a s ) GeV, and u(a s ) « (0.35 + a s ■ 1.2) GeV. This 
behaviour is typical for other hadrons, see [T1J for the discussion of glueballs 
and II for a review of hadronic spectra. 

^From dimensional point of view the spectrum corresponds to the increas- 
ing dimension of valence gluon operators a iy diCik, didkdi etc. (and not to the 
dimension of operators E^B^ DiB k etc.). 

Note that the lowest level, 1 , is of the electric type and there is a gap 
of 0.5 -T- 0.6 GeV between the ground and excited level, while the distance 
between L = and L = 2 levels is equal to 0.73 GeV (a s = 0) or 1 GeV 
(a, = 0.15). 

It is instructive to compare our spectrum to the gluelump spectrum of 



the MIT bag model ||, [24]. 



For A 1 / 4 = 0.315 GeV and a, = 0.23 the masses obtained in P41 



are 



M(l+-) = 1.43GeV,M(2~ ) = 1.97GeV, M(l~ ) = 1.98GeV 

M(3 + ~) = 2.44GeV, M(2 + ") = 2.64GeV. (26) 

One can see the inverse ordering of the first levels in (E^), namely magnetic 
level 1 H in MIT bag model is below the electric one 1 and the gap between 
them is approximately the same as in QCD string model, but the ordering 
is reversed. 



9 



Here the difference between two models is most pronounced (in contrast 
the glueball spectrum, where the lowest two-gluon states have the same or- 
dering in both models, see JIT] and [2"5[|). Hence the independent check of 
lattice calculations can in principle decide which of the models is closer to 
reality. 

Now the recent lattice calculations |2f| (see also for earlier lattice 
studies) yield the mass levels shown in Table 5. One should note at this 
point, that the (divergent) mass renormalization terms are not subtracted 
from the data and therefore absolute scale is missing. We therefore have 
placed for the sake of comparison in Table 5 the lowest l + ~ level at the same 
mass as given by our calculations, see Table 2 and 4 for 6t s = 0.15. 

Another feature of data from is that only spacial links have been 
used for gluelump operators at initial and final states, and hence the states 
in Table 4 containing Ej are not directly excited in measurements in |2"6" 



(however these states could in principle be excited indirectly due to mixing, 
i.e. with smaller overlap). The states in Table 4 which are excited directly 
in |26] are marked by asterix and compared with lattice results in Table 5. 

One can see in Table 5 a good correspondence between the levels, and an 
approximate mass degeneration of 2 + ~,3 +_ levels is present in both results. 

However the 1 ,2 levels, almost degenerate in the QCDSM have a 
gap of 0.22 GeV on the lattice. This might be explained by the mixing of the 

1 level with the ground state, having the same quantum numbers, whereas 

2 has no low- lying counterpart. 

As was mentioned above, the crucial test of the QCDSM is the search for 
the ground state level 1 , which should lie 0.5 -r- 0.6 GeV below the lowest 
l + ~ level found on the lattice. If this level is not found, it would possibly 
mean that gluelumps are more like bags than strings. However no deriva- 
tion of the bag model from QCD was ever made and theoretical grounds for 
this model are of intuitive character, whereas the QCDSM follows directly 
(with few assumptions supported by lattice data) from the QCD lagrangian. 
Therefore the existence of the 1 ground state for the gluelump seems very 
plausible. In conclusion of this section we discuss shortly an important con- 
nection between gluelumps and field strength correlators D(x), Di(x) intro- 
duced in JT4J and measured on the lattice in [[H]] using cooling procedure 
and in [27[] using RG smoothing. It is clear from Table 4 that two lowest 
states, 1 and l + ~, correspond to the electric and magnetic correlators re- 
spectively, and therefore the masses M(l ) = 1.12 GeV and M(1 H ) = 1.7 
GeV (a s = 0.15) should give the exponential slopes of those correlators (the 
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notation Dm and D^ was also used in [17 coiciding with electric and magnetic 
respectively). The lowest mass of M « 1 GeV is in agreement with mea- 
surements in ||T7J , however for D± the same slope was found in contrast with 
present calculations. There are no general symmetry arguments, why masses 
(slopes) in both functions should be the same, (in contrast to electric and 
magnetic condensates at zero temperature) and this disagreement between 
gluelump masses and slopes in D^,D± calls for further investigation. 

One should note in addition, that the independent calculation of the 



slope from the selfcoupled equations for correlators made in | 28fl , also gives 



the lowest (electric) mass around 1 GeV, in agreement with the present cal- 
culations for (a s = 0.15). An analysis for the same quantity in the framework 
of the QCD sum rules |2Pj gives a larger value in the quenched case, while for 
the magnetic mass the result was unstable. Finallyin the lattice study |3(| 
of electric and magnetic field correlators made without cooling or smooth- 
ing procedures,electric and magnetic masses are obtained in good agreement 
with present calculations, specifically magnetic mass turned out to be around 
0.5 GeV heavier. 



6 Conclusions 

Results of the present calculations in the QCDSM framework are in general 
agreement with recent lattice data [[^J , using the restricted set of gluonic 
operators, and with previous lattice and analytic calculations of the slope of 
the Dm function. It is not clear from the gluelump spectrum point of view 
why the correlator D± (magnetic field correlator) should have the same slope 
as that of Dy, as it was found on the lattice |17| . With the lowest gluelump 
mass equal to 1.2 -=-1.5 GeV in the present calculations, the onset of the 
adjoint string breaking appears at the string length ro = 1.2 4- 1.5 fm, in 
agreement with recent lattice data ||18|| . 

The fundamental role of gluelump states in the theory of the QCD vacuum 
makes it important to resolve the existing problem of the spectrum and of 
the relation with vacuum correlators D±, Dy. 
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Table 1 



Valence gluon masses no(n,L) (upper entries) and reduced eigenvalues 
a(n,l) (lower entries). a f = 0.18 GeV 2 ,a adj = \<j } = 0.405 GeV 



L\n 





1 


2 





0.746 


1.135 


1.422 




2.3381 


4.0879 


5.520 


1 


0.98 


1.297 


1.554 




3.3613 


4.8845 


6.216 


2 


1.168 
4.248 


1.443 
5.63 




3 


1.33 
5.053 







Table 2 

Unperturbed masses M (n,L), string correction AM L and spin averaged 
masses M(n, L) for L — 0, 1, 2, 3 and n = (upper entries and two last lines) 
and n = 1 (lower entries). 



\L 





1 


2 


3 


M 




1.492 


1.96 


2.336 


2.66 


(GeV) 




2.27 


2.6 


2.886 




AM L 







-0.0581 


-0.103 


-0.139 


(GeV) 







-0.025 


-0.0546 




M 




1.492 


1.9 


2.233 


2.52 


(GeV) 




2.27 


2.575 


2.83 




M(a s = 


0.15) 


1.12 


1.73 


2.15 


2.46 


M(a s = 


0.195) 


0.982 


1.67 


2.10 
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Table 3 

Perturbative and Thomas spin-orbit martix elements defined as in ( |TT| - |i"3|) 



L 


1 


2 


3 


foipert) 


0.033 


0.0096 


0.00425 


^(Thomas) 


-0.0025 


-0.0024 


-0.0021 



Table 4 

Gluelump masses for lowest states J PC are listed together with gluonic 
state operators in the background field formalism and in the general nota- 
tions. Asterix marks the states which are directly generated on the lattice 
PI 



State 


L 


n 


operator 


Mass 
(GeV) 


backgr. 


general 


\ — 








Oi 


Ei 


1.492 


1+- 


1 







B t 


1.87* 


2+- 


1 





\diQ>k) symm 


(.Dj E k ) symm 


1.93 


\ — 





1 


Oi 


Ei 


2.27 


\— 


2 





{did k ai)j 


DiB k e ik i 


2.21* 


2~ 


2 





n 


\D% B k J S ymm 


2.226* 


3~ 


2 





n 




2.24 


2+- 


3 





(did k dia m )j 


{DiD h Bi)j 


2.51* 


3+- 


3 





n 


» 


2.52* 


4 +- 


3 





>! 


(DiDkDiEm)^ 


2.53 


1+- 


1 


1 


{dia k )i 


B t 


2.57* 


2+- 


1 


1 


{dia k ) 2 




2.57 


0++2++ 








K a i a k )J 




2.61* 
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Table 5 



The gluelump states generated by spacial operators, with masses com- 
puted analytically in the present paper (a s = 0.15) and masses computed on 
the lattice [p6|] . For better comparison results of Table 2 corrected for larger 



(jj = 0.22 used on lattice. 



jPC 


1+- 


1~ 


2— 


2+- 


3+- 


0++ 


M (GeV) 
lattice 


1.87 


2.23 


2.45 


2.84 


2.84 


2.96 


M (GeV) 
QCD string 


1.87 


2.34 


2.36 


2.70 


2.71 


2.78 
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Appendix 



Two-gluon gluelumps 

Consider a static source at point r3 = and two gluons at r = ri and 
r = r 2 , connected by a fundamental string, passing through these three 
points. The Hamiltonian for the case, when both orbital momenta are zero 
(no string rotation correction) can be written in analogy to (0) as 

2 2 

hQii,to) = + + + *-,{|r a | + |r 2 | + |n - r 2 |} (A.l) 



This 3 body problem can be treated by the hyperspherical formalism pi 
Introducing the hyperradius p and denoting pi = /i 2 = p, one has 

P 2 = rf + d\ r) = ri2 /V2, ^ = ( ri + r 2 )i= (A.2) 

one obtains the Hamiltonian for the given grand orbital momentum K = 
0,1,2,... 

h = + » - + + UI<{P) (A ' 3) 



where Ux(p) for K = is the potential term in (|A.l ) averaged over all 



hyperangles. Using standard expressions from [31| one obtains 



tt(\ r r 32y^(l + y/2) 

U {p) = C ap, C = — (A.4) 

157T 



Eq. ( |A.3|) can be solved numerically and results are tabulated, but it appears 



that the accuracy of about 1% for the eigenvalue can be obtained by the 
stationary point method applied to the effective potential W(p) 



W(p)= {K + f ( * +V +U K ( P ); W'(p = p ) = 0; (A.5) 
2pp 2 



while the second derivative at p = p defines the radial excitation energy uj: 

w(J^p£ s ^(p-p^ (A . 6) 
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In this way one obtains the mass eigenvalue M(/x), and finally the stationary 
point condition for /x, = /xo) = defines gluon constituent mass /x 

and the resulting mass eigenvalue M(/xo), given in (p2j)-(|24]). One can note 
that radial excitations are given by ujq = 0.955 GeV, and "grand-orbital 
excitation" from K = to K = 1 is given by the interval AM = 0.712 GeV. 
The hyperfine splitting is given by the Hamiltonian (cf the discussion in 



TTf for 3g glueballs) 



H 88 = 3W ^MT* ia ' (S® (A. 7) 



The last factor in ( |A.7|) is easily computed 



(5 (3) (ri 2 ))p = ^|, Po = 1.15/aiq. (A.8) 



This yields the spin splitting of the levels 



AM SS = S (1) S (2) 0.49/x ^cv (A.9) 



For a, = a, = 0.15 one obtains 



AM SS = 0.0525GeV ( j = ° 2 ) ■ (A.10) 



In this way one gets the mass M(0 ++ ) in Table 5. 
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